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(57) Abstract 

The aim of the invention is, in particular, to control and/or 
regulate a recovery boiler used for the production of cellulose, 
whereby signals are determined from the process data produced 
for the input and output signals. More specifically, to this end raw 
data are entered into a data preprocessing unit, the processed data 
are used for system identification, and the desired target values 
for output signals, particularly for set values, are inputted by an 
operator. An optimization strategy is derived, which proposes 
a change in the input signals. To this end, the structured data 
are introduced into the regulating process. The associated device 
comprises an additional device (advisory system) for generating 
signals. In this way, the operator is presented with a prediction 
which corresponds to the control proposition and process model. 
The process model is re— corrected, i.e. adapted, in accordance 
with current process data. 

(57) Zusammenfassung 

Insbesondere soli ein Ruckgewinnungskessel (sog. Re- 
covery-Boiler) bei der Zellstoffherstellung gesteuert und/oder 
geregelt werden, wobei aus den flir Eingange und Ausgange an- 
fallenden ProzeBdaten Sign ale ermittelt werden. Im einzelnen 
werden dazu die Rohdaten in eine Einheit zur Datenvorverar- 

beitung gegeben, die aufbereiteten Daten fur eine Systemidentinzierung verwendet, und vom Operator die gewunschten Zielwerte der 
Ausgange, insbes. Sollwerte, eingegeben. Es wird eine Optimierungsstrategie abgeleitet, die ein Veranderung der Eingange vorschlagt, 
wozu die strukturierten Daten in den RegelprozeB eingegeben werden. Bei der zugehdrigen Vorrichtung ist eine zusatzliche Einrichtung 
(Advisory-System) zur Generierung von Signalen vorhanden. Damit wird dem Operator entsprechend dem Steuervorschlag und dem ProzeB- 
modell eine Pradiktion angeboten. Das ProzeBmodell wird aufgrund von aktuellen ProzeBdaten im Sinne einer Adaption nachkorrigiert 




LEDIGLICH ZUR INFORMATION 



Codes zur Identifizierung von PCT-Vertragsstaaten auf den Kopfbogen der Schriften, die Internationale Anmeldungen gemass dem 
PCT veroffentlichen. 



AL 


A 1 ban i en 


ES 


Spanien 


LS 


Lesotho 


SI 


Slowenien 


AM 


Armcnien 


FI 


Finnland 


LT 


Litauen 


SK 


Slowakei 


AT 


Osterreich 


FR 


Frankreich 


LU 


Luxemburg 


SN 


Senegal 


AU 


Australien 


GA 


Gabun 


LV 


Lett land 


SZ 


Swasiland 


AZ 


Aserbaidschan 


GB 


Vereinigtes Kdnigreich 


MC 


Monaco 


TD 


Tschad 


BA 


Bosnien-Herzegowina 


GE 


Georgien 


MD 


Republik Moldau 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagaskar 


TJ 


Tadschikistan 


BE 


Belgien 


GN 


Guinea 


MK 


Die ehemalige jugoslawische 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Griechenland 




Republik Mazedonien 


TR 


TQrkei 


BG 


Bulgarien 


HU 


Ungam 


ML 


Mali 


TT 


Trinidad und Tobago 


BJ 


Benin 


IE 


Irland 


MN 


Mongolei 


UA 


Ukraine 


BR 


Bras i lien 


IL 


Israel 


MR 


Mauretanien 


VG 


Uganda 


BY 


Belarus 


IS 


Island 


MW 


Malawi 


US 


Vcreinigte Staaten von 


CA 


Kanada 


IT 


Itatien 


MX 


Mexiko 




Amerika 


CF 


Zentralafrikanische Republik 


JP 


Japan 


NE 


Niger 


uz 


Usbekistan 


CG 


Kongo 


KE 


Kenia 


NL 


Niederlande 


VN 


Vietnam 


CH 


Schweiz 


KG 


Kirgisistan 


NO 


Norwegen 


YU 


Jugoslawien 


CI 


Cfite d'lvoire 


KP 


Democrat is c he Volksrepublik 


NZ 


Neuseeland 


zw 


Zimbabwe 


CM 


Kamerun 




Korea 


PL 


Polen 






CN 


China 


KR 


Republik Korea 


FT 


Portugal 






CU 


Kuba 


KZ 


Kasachstan 


RO 


Rum an i en 






CZ 


Tschechische Republik 


LC 


St. Lucia 


RU 


Russische FO deration 






DE 


Deutschland 


LI 


Liechtenstein 


SD 


Sudan 






DK 


Dan cm ark 


LK 


Sri Lanka 


SE 


Schweden 






EE 


Estland 


LR 


Liberia 


SG 


Singapur 







WO 98/27474 



PCT/EP97/07096 



1 

Beschreibung 

Verfahren und Vorrichtung zur Prozefifuhrung einer verfahrens- 
technischen Anlage 

Die Erfindung bezieht sich auf ein Verfahren zur Prozefifuh- 
rung einer verf ahrenstechnischen Anlage, insbesondere eines 
Riickgewinnungskessels (sog. Recovery-Boiler) bei der Zell- 
stof fherstellung, wobei aus fur Eingange und Ausgange an- 
fallenden Prozefidaten Signale ermittelt und zur Steuerung 
und/oder Regelung des Prozesses verwendet werden. Daneben 
bezieht sich die Erfindung auch auf die zugehorige Vorrich- 
tung, die vorzugsweise zur Steuerung und/oder Regelung des 
Recovery-Boilers als verf ahrenstechnische Anlage geeignet 
ist, wobei der Anlage wenigstens eine Steuer- und/oder 
Regeleinrichtung zugeordnet ist. 

Der Recovery-Boiler in einer Zellstof f f abrik hat im wesent- 
lichen zwei Aufgaben: 

1. Riickgewinnung von Chemikalien aus dem anorganischen Anteil 
der Schwarzlauge 

2. Produktion von Dampf aus dem organischen Anteil der 
Schwarzlauge zur Energiegewinnung . 

Der HerstellungsprozeS fur Zellstoff beinhaltet das Kochen 
von Holzschnitzeln, wobei nach dem Kochungszyklus ein Zell- 
stoffbrei (Pulpe) und eine schwache Schwarzlauge (weak black 
liquor) entstehen. Die schwache Schwarzlauge (ca. 12 bis 18 % 
feste Stoffe) wird durch Eindampfung zu einer starken 
Schwarzlauge (strong black liquor = BL) mit 50 bis 55 % 
Feststoffen konzentriert . Durch einen weiteren Verfahrens- 
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schritt wird die Schwarzlauge (BL) auf 60 bis 80 % Feststoffe 
konzentriert, urn die Verbrennung zu erleichtern. 

Letztere Schwarzlauge (BL) wird in den Recovery-Boiler ein- 
gespruht, wobei BL-Tropf chen entstehen. Diese BL-Tropf chen 
verbrennen im Recovery-Boiler und produzieren u.a. C0 2/ Na 2 S0 4 
und ggfs. auch Kohle. Die organischen Bestandteile verbrennen 
zum groSen Teil im Gasraum; eine weitere Umsetzung erfolgt in 
der Kohleschicht auf der Oberflache eines Schmelzbettes 
(char-bed) , das sich am Boden des Recovery-Boilers ausbildet. 
Durch die anorganischen Bestandteile der Schwarzlauge ent- 
steht eine Sodaschmelze, die aus dem Recovery-Boiler ab- 
flieSt. Diese Sodaschmelze wird in nachf olgenden Schritten 
des Riickgewinnungszyklus (Recovery cycle) weiterverarbeitet 
und dem KochprozeS wieder zugef uhrt . 

Das sog. RuSblasen (Sootblowing) ist ein notwendiger 
Reinigungsvorgang der Warmeuber gangs f lachen im Recovery- 
Boiler. Dabei wird ein gewichtiger Anteil des im Recovery- 
Boiler erzeugten Dampfes verbraucht (ca. 5 bis 10 %) . Das 
Rufiblasen wird zyklisch durchgef uhrt . Daher mu£ u.a. zwischen 
der Ef f izienzsteigerung des Recovery-Boilers durch das 
Sootblowing und dem dafur notwendigen Dampf verbrauch 
optimiert werden. 

Aufgabe der Erfindung ist es, ein Verfahren und eine Vor- 
richtung anzugeben, mit denen vorbeschriebener Prozefi in 
geeigneter Weise gefuhrt werden kann. 

Zur Losung dieser Aufgabe wird erf indungsgemaS das im 
Patentanspruch 1 angegebene Verfahren insbesondere zwecks 
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Optimierung der ProzeSf uhrung in der Zellstoff- und Papier- 
industrie vorgeschlagen . Vorteilhafte Weiterbildungen sind in 
Unteranspriichen angegeben, wobei insbesondere die Anspruche 9 
bis 11 praxisrelevante Ausbildungen fur einen Recovery-Boiler 
5 kennzeichnen. Die zugehorige Vorrichtung wird in den Patent - 
anspruchen 12 und 13 angegeben. Eine solche Vorrichtung ist 
geeignet, dem Anlagenf ahrer als Operator eines anlagen- 
technischen Systems geeignete Hilf estellung zu geben. 

10 Bei der Erfindung ist besonders vorteilhaft, daS uber das 

Advisory -System dem Operator konkretere Steuervorschlage zum 
Fahren der Anlage gemacht werden konnen und dariiber hinaus 
sogar eine Pradiktion als Vorhersage uber den Anlagenbetrieb 
moglich ist. Dabei kann das dazu notwendige Prozefimodell 

15 aufgrund von aktuell en ProzeSdaten nachkorrigiert werden, so 
dafi das Modell dem aktuellen Anlagen-Zustand entspricht. 

Weitere Einzelheiten und Vorteile des erf indungsgemaiSen Ver- 
fahrens und der zugehorigen Vorrichtung speziell zur ProzeiS- 
2 0 f uhrung eines Recovery-Boilers ergeben sich aus nachf olgender 
Figurenbeschreibung . Es zeigen 

Figur 1 eine schematische Darstellung eines Recovery-Boilers, 

Figur 2 die Grobstruktur des zugehorigen Regelkreises , 
25 Figur 3 den Aufbau des Advisory-Systems und 

Figur 4 ein Kamerabild eines Schmelzbettes des Recovery- 
Boilers gemaS Figur 1 und 

Figur 5 die Bedienoberf lache eines Interface des Advisory- 
Systems fur die Mensch/Maschine-Kommunikation. 
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Die Figuren werden nachfolgend teilweise gemeinsam beschrie- 
ben. Die Erfindung wird daran beispielhaft speziell fur einen 
Recovery-Boiler erlautert. Sie ist aber auch fur andere ver- 
f ahrenstechnische Anlagen insbesondere in der Zellstoff- und 
Papierindustrie einsetzbar. 

Der eingangs bereits erlauterte Verbrennungs- und 
RuckgewinnungsprozeS bei der Zellstoff herstellung kann durch 
die maSgeblichen vier Eingange (Prozefiimpuls = Manipulated 
variables) u beeinflufit werden: 

1- Menge der einstromenden Schwarzlauge (BL-Flow) 

2. Zugefuhrte Primarluft (Primary Air = PA). 

3 . Zugefuhrte Sekundarluf t (Secondary Air = SA) . 

4. Zugefuhrte Tertiarlauft (Tertiary Air = SA) . 

5. Druck der einstromenden Schwarzlauge (BL-Pressure) . 

Die wesentlichen Ausgange (ProzeSoutputs) y lassen sich 
speziell bei einem Riickgewinnungskessel in der Papier- und 
Zellstoff industrie durch folgende Variablen beschreiben: 

1. Schmelzbett (char bed.) 

Char-bed - Gestalt und zwar Hohe und Form 

Char-bed - content, d.h. der chemische Gehalt der Einzel- 

komponenten 
Char-bed - Temperatur 

Der Bettgestalt kommt besondere Bedeutung zu. 

2 . Gas : 

NOx - Gehalt 

SOx - Gehalt 

CO - Gehalt des Gases 

02 - Uberschufi (Excess 02) im Gas 
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H2S + xx - Total reduced Sulfur (TRS) 

3. Tropf chengroEe (Droplet size) 

4 . Zahl und Volumen der zu den Superheatern vom Gasstrom 
hinauf gerissenen Partikel (Carry Over) 

5 . Dampf : 

Dampf durchf lulSmenge F 
Dampfdruck P 
Dampf temperatur T. 

6. Temperatur im Kessel (Temperaturprof il) 

7. Reduktionsgrad der Grunlauge 

8 . pH-Wert der Asche 

Bild 1 zeigt ein vereinf achtes Schema eines Recovery Boilers 
mit seinen Ein- und Ausgangen. Dargestellt ist ein Kessel 1 
mit Einlassen 2 bis 6 sowie Auslasse 7 und 9. Im einzelnen 
wird uber die Einlasse 2 und 3 Schwarzlauge (BL) eingespruht, 
aus der sich nach Verbrennung eine Schmelze in einem Schmelz- 
bett 10 variabler Gestalt bildet. Verbrennungsluf t wird liber 
die Einlasse 4 bis 6 zugefiihrt. Uber den AuslaS 7 konnen 
Proben der Schmelze zur Reduktionsgradbestimmung entnommen 
we r den. 

Im oberen Bereich 8 des Kessels 1 wird mit der Verbrennung 
Dampf erzeugt und uber einen GasauslaS 9 ausgegeben. Es sind 
MeSstellen fur die Dampf variablen Durchflufi (F) , Druck (P) 
und Temperatur (T) vorhanden. Weiterhin werden die Gas- 
bestandteile wie 0 2/ CO, S0 2 , H 2 S und NO x bestimmt . Auch der 
pH-Wert der anfallenden Asche kann erfaSt werden. 

Die allgemeinen Ziele der ProzeSsteuerung beim Betrieb des 
Recovery-Boilers gemaS Figrur 1 sind 
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1. Verbesserung der Ruckgewinnung der anorganischen 
Chemikalien (chemical recovery) 

2 . Erhohung der Energiegewinnung aus organischen Anteilen 
(energy conversion efficiency) 

3. Erhohung der Verf ugbarkeit des Recovery Boilers (Betriebs- 
kosten) 

4. Einhaltung der gesetzlich festgelegten Maximal -Konzen- 
trationen von Schadstoffen in den Emissionen. 

Letztere Ziele konkurrieren xiblicherweise miteinander, so daS 
bei der Steuerung eines Recovery-Boilers eine komplexe Pro- 
blematik vorliegt. Gegenwartig werden Recovery-Boiler vom 
Anlagenf ahrer , der nachfolgend als Operator OP bezeichnet 
wird, iiberwacht und gesteuert . Der Operator OP fahrt den 
ProzeS entsprechend seiner Erfahrung und insbesondere anhand 
von vorhandenen ProzeSrichtlinien . 

Zur Einstellung der StellgroEen u= (PA, SA, TA, BL-Flow, BL- 
Press) gibt es eine Grundaut omat i s ierung , die ublicherweise 
mit sog. Low-Level - Controllern durchgefiihrt wird , was hier 
als gegeben vorausgesetzt wird. Uber der Ebene der Low-Level - 
Controller ist nunmehr eine weitere Regelungs-/Optimierungs- 
ebene vorgesehen, die 

1. voneinander unabhangige Ausgange (ProzeSoutputs) y u auf 
vorgegebene Sollwerte (Setpoints) y d einstellt und 

2. voneinander abhangige Ausgange (Prozefioutputs) p a mit 
Hilfe einer entsprechenden Gewichtung Q so bewertet, daf£ 
sich im stationaren Fall gewunschte Verhaltnisse zwischen 
den Komponenten von p o einstellen. 



WO 98/27474 



PCT/EP97/07096 



Eine entsprechende Grobstruktur des Regelkreises zeigt Figur 
2. In einer Einheit 20 wird ein Prozefi 21, der beispielhaft 
durch den Betrieb des Recovery-Boiler aus Figur 1 definiert 
ist, durch einen Controller 22 geregelt. ProzeSdaten werden 
5 auf ein sog. Advisory- System 3 0 vibergeben, das dem Operator 
OP Hinweise zum Einrichten der Setpoints sp am Controller 22 
liefert. 

Fur die Reglerkonzeption und die zugehorige Optimierung ist 
10 eine hinreichende ProzeSkenntnis notwendig, die z.B. auf eine 
datengetriebene Identif ikation eines Prozefimodells auf baut . 
Das entsprechende ProzeSmodell wird dabei aus Betriebsdaten 
gewonnnen. Im folgenden wird zunachst das ProzeJSmodell 
beschrieben. Danach folgt eine Beschreibung der Schritte 
15 a) Datenvorverarbeitung (data reconcilliation) 

b) Systemidentif ikation (system identification) 

c) Optimierung 

d) Operatorberatung (operator advice) . 

20 Aus Figur 3 ist der Verf ahrensablauf im Advisory- System 30 
aus Figur 2 naher ersichtlich. Dazu sind Ablauf einheiten 32 
bis 3 6 zusammengeschaltet , die den sof twaremaSigen Verf ah- 
rensablauf wiedergeben. Die Rohdaten. gelangen uber einen 
Eingang 31 auf die Einheit 3 2 zur Datenauf bereitung und von 

25 dort auf die Einheit 33 zur Systemidentif ikation . Hier 

erfolgt gleichermaSen eine Daten-Clusterung . Der Einheit 34 
zur Optimierung ist weiterhin ein externer Optimierer 3 5 
zugeordnet . Von dort werden die Daten auf die Einheit 36 zum 
Control Advice fur den Operator OP aus Figur 2 gegeben. 
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Zur Prozefiidentif ikation und Regelstrategie wird zunachst das 
ProzeSmodell betrachtet : Das dynamische ProzeSmodell hat im 
diskreten Fall folgende allgemeine Form 

x k+ i = f(x k ,u k ) (la) 

y k = i(x k ,u k ) db) 

mit 

x k - Zustandsvektor 

u k - Stellgrofienvektor 

y k - Outputvektor 

Im vorliegenden Fall werden folgende Annahmen getroffen: 

1 . Die Prozefioutputs und ZustandsgroEen werden gleichgesetzt 

Yk ~ X k 

2. Das System wird als nichtlineares , aber statisches 
Ubertragungsglied angesehen 

y k =i(S k ). (2) 

Diese Annahme liber die statische Ubertragungseigenschaf t ist 
darin begriindet, daS 

1. das Open-Loop-System gemaS Gl'en (la) und (lb) als stabil 
anzunehmen ist und 

2. die groSte Systemzeitkonstante als klein gegenuber der 
Abtastzeit angenommen wird. 

Letzteres bedeutet, daS im Augenblick einer Messung bzw. 
eines Stelleingrif f s sich das System immer im eingeschwun- 
genen Zustand bef indet . Dieses soil aus Grunden der Verein- 
fachung der Regelungsstrategie naherungsweise auch fur die 
Bettgestalt gelten, worauf weiter unten eingegangen wird. 
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Eine Besonderheit des Systems ist, daS es weniger Eingange 
als Ausgange hat : 

dim(u) < dim(y) . (3) 

Das bedeutet im vorliegenden Fall, daS mit den vorhandenen 
Eingangen die Ausgange nicht beliebig beeinfluSt werden 
konnen, was fur die Systemidentif ikation zwar keine Rolle 
spielt, aber fur die Regelstrategie von Belang ist. 

Im folgenden wird das System gemaJS Gl . (2) um einen 
beliebigen Arbeitspunkt u k linearisiert 

A y k = J c'Au k , (4) 



mit Ay k = y k+1 -y und der Jacobi -Matrix J = — 



Gleichung 



(4) beschreibt somit einen linearen Zusammenhang zwischen 
AnderungenAu k am Eingang des Systems und den daraus 
folgenden Anderungen Ay k am Ausgang des Systems. 
Fur ein Teilsystems des Recovery-Boilers gemaS Figur 1 seien 
die Eingangsanderung Au k =ABL k und die dazugehorenden Aus- 
gangsanderungen Ay k = (A0 2 , ACO) k T . Die entsprechende Jacobi- 
Matrix (Jacobian) lautet 



J = ( 



<?B L 



dC 0 



B L = B L 1 L 



y 

\ L = B L k / 



< 5 > 

BL = 



Die Identif ikationsauf gabe besteht nun darin, fur bestimmte 
Arbeitspunkte u 0 die Jacobi-Matrizen zu ermitteln, wahrend 
die Regelungsaufgabe darin besteht, fur die ermittelten 
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Jacobi-Matrizen eine Regelstrategie zu entwerfen, die das 
System in den Arbeit spunk ten u G optimiert. 

Zur Datenvorverarbeitung (Data reconciliation) wird bei den 
anfallenden Daten vom stationaren ProzeS ausgegangen, so daS 
die Daten von instationaren Vorgangen entfernt werden konnen. 
Die ProzeSdaten konnen stark verrauscht sein, insbesondere 
was die Gasvariablen betrifft. Die Datenvorverarbeitung be- 
handelt daher verrauschte Prozefidaten, die aufierdem aufgrund 
von Storungen in der Sensor ik sehr stark aus einem vorge- 
gebenen Toleranzband herausf alien konnen. 

Zur Systemidentif ikation kann zunachst eine datengetriebene 
Systemidentif ikation mit Clusterung erfolgen. Das Ziel einer 
Clusterung ist, aus vorliegenden Daten die Arbeitspunkte 
herauszuf inden , 

Die Arbeitspunkte werden im niedrigdimensionalen Eingaberaum, 
also insbesondere fur die Eingange u = (PA, SA, TA, BL- 
Flow) , definiert. Diese Besonderheit des Verfahrens wurde 
deshalb gewahlt, da wegen der geringeren Anzahl der Eingange 
entsprechend Gl . (3) einige Ausgange voneinander abhangig sind 
und daher eine unabhangige Clusterung des Ausgaberaums nicht 
moglich ist . 

Die Daten-Clusterung erfolgt mittels eines bekannten Fuzzy- 
Clustering-Algorithmus , der in folgenden Schritten arbeitet : 
1. Vorgabe der Merkmale, nach denen geclustert wird. Im 

vorliegenden Fall sind dieses die Eingangsvariablen u = 
(PA, SA, TA, BL-Flow) 
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2. Vorgabe der Anzahl n der zu erwartenden Cluster. Bei den 
vorliegenden Daten hat sich eine Clusteranzahl von n=10 
als am besten erwiesen. 

3 . Berechnung der Clusterzentren 

5 4. Ermittlung der Zugehorigkeitsf unktionen fur u = (PA, SA, 
TA, BL-Flow) , die fur jede Komponente von u angeben, 
welchen Zugehorigkeitswert zu einem Cluster Cjein Datum 
besitzt, das einen bestimmten Abstand d von dessen 

Clusterzentrum uj. hat. Diese Zugehorigkeitswerte werden 
10 bei der Systemidentif ikation als Datengewichte 

beriicksichtigt . 

Anschliefiend erfolgt eine Identif ikation des Prozefimodells 
auf der Grundlage derjenigen Daten, die nach der Datenvor- 
15 verarbeitung vorliegen. 

Die Identif ikation des Prozefimodells dient der spateren 
Vorhersage (Pradiktion) der Ausgange (Prozefioutputs) im 
Advice -Teil der Steuerung fur den Recovery-Boiler. Bei der 

20 Modellbildung wird von einem statischen Prozefi ausgegangen, 
bei dem jede innere Prozefidynamik vernachlassigt wird. Das 
ist deshalb gerechtf ertigt , da die Zeitkonstanten schneller 
gasdynamischer Vorgange aber auch langsamer thearmischer 
Prozesse unterhalb der Eingrif f szeit der Operators liegen. 

25 Die Modellierung geht daher zunachst von einem nichtlinearen, 
aber statischen Zusammenhang zwischen Inputs wund Outputs y 
des Prozesses aus : 

y = fm(M)- (6) 
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Diese Funktion wird in verschiedenen Arbeit spunk ten des 
Input -Raumes linearisiert und mit Hilfe des bekannten 
Takagi/Sugeno-Modellansatzes approximiert . Es ergibt sich 



c 



y = ^M t (") <Jj ■ u + const 1 ) 



(7) 




mit const' = — J* • u' + y* 

(w) - Zugehorigkeitsgrad (degree of membership) of des i. clusters 
JJ - Jacobian des i. cluster 



const 1 



Konstante des i. clusters 



Kern des i. clusters 
Ausgangsvektor zu u 1 



Die Identif ikation benutzt ein rein datengetriebenes 
Verfahren, bei dem zunachst eine Clusterung des Input -Raumes 
vorgenommen wird und anschliefiend die Jacobians sowie die 
10 konstanten Biases identif iziert werden. 

Die dazu benutzten Datensatze basieren auf Tagessatzen von 
ProzeSdaten aus drei bis vier Monaten, die folgender Prozedur 
unterzogen werden : 

15 

a) lineare Interpolation der sog. Infrequent Measurements 

b) Einfiigen von extra gemessenen Bettdaten 

c) Herausstreichen von Datensatzen, deren Inputs oder Outputs 



20 



e) 



d) 



maximale oder minimale Grenzen iiber- bzw. unterschreiten 
Filtern der Outputs Y_CUR zu Y_FILT 

Ausschneiden von Transienten (15 min nach einer Input - 
Anderung der Luftstrome und maximal 10 mal in einem 
Tagessatz) 
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f ) Zusammenf assen der Tagesdatensatze zu einem groBen 
Datensatz 

g) Gegebenenf alls lineare Interpolation der Bettdaten 

h) Komprimierung der groSen Datensatzes mit einem frei 
5 wahlbaren Kompressionsf aktor 

Der Algorithmus zur Clusterung ist folgendermafien struktu- 
riert : Ausgehend von einer vorgegebenen Clusteranzahl c und 
einer Unscharfe m werden die Clusterzentren bis zu einer 
10 maxitnalen Iterationsanzahl maxiter bzw. bis zu einer Ab- 

bruchschwelle u_eps iterativ berechnet . Mit den Zugehorig- 
keitsgraden u und den Clusterzentren v. ergeben sich dann 
f olgende Schritte : 



15 



1. Berechnung der Clusterzentren v, 



20 



v =-^- 



-, Vi e{l, — ,c} , 



(8) 



2>* 



2 . Berechnung einer Fuzzy-Scatter-Matrix S t 
S,-=2>* •(x A -v,).(x jk -v / ) r ,Vi€{l t ...,c} , 



(9) 



Jfcal 



3 . Extrahieren der Eigenvektoren d t der maximalen 

Eigenwerte jeder S t , 
4 . Berechnung der Zugehorigkeitsgrade u ik und Abstande 

D ik zwischen den Punkten und den Clusterzentren v, : 



A* = - v,|| 2 - a • 5J - v,) r • d 9 ) und u ik = 



>=1 



m-l 



(10) 
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Vie{l,-,c>, Vke{l,.-,N} 
bis {l=maxi ter oder 

Die anschliefiende Regression berechnet fiir jedes Cluster die 
Jacobian und den konstanten Bias anhand der Fuzzy- Clustering - 
Ergebnisse und der gereinigten und normalisierten Eingangs- 
variablen sowie jeder einzelnen gereinigten und normalisier- 
ten Ausgangsvariablen. In diesem Zusammenhang bedeutet 
„gereinigt u , daS sogenannte AusreiSer in den Inputs um einen 
Medianwert herum entfernt werden, wobei in diesem Fall sowohl 
der Input- als auch der Output -Vektor gestrichen werden. 

Vor jeder Regressionsberechnung fur ein bestimmtes Cluster 
wird jeder Datenpunkt durch seinen Zugehorigkeitsgrad zu 
diesem Cluster gewichtet . 

Das Modell wird derart berechnet, daS jeder einzelne Ausgang 
getrennt berucksichtigt wird. Auf diese Weise erhalt man ein 
sog. MISO-Modell. Der Vorteil dieser Vorgehensweise ist, daS 
fur die Erstellung eines Modells fur einen speziellen Ausgang 
die maximale Anzahl der zur Verfugung stehenden Daten benutzt 
wird. 

Der Algorithmus ist f olgendermaSen strukturiert : 
1 ♦ Es seien 



14 

max(||4 - u'- l \\) <u_eps} . (n) 



[4 = MiklMik) t>zw. 



(12a) 
(12b) 
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die gewichtete Eingangsmatrix bzw. der gewichtete Ausgangs- 

vektor fur das Cluster k. [X]J 1st eine n*n_in Matrix und 

ist ein n*I Vektor, wobei n_in die Dimension des Eingaberaums 
und n die Anzahl der Datenpunkte sind. Der normalisierte 
Regressionsvektor b jk wird dann f olgendermaSen def iniert : 

wobei der letzte Koefficient von b Jk den konstanten Koeffi- 
zienten der Regression bezeichnet (Bias) . Fur jeden Ausgang 
ist die Zahl n unterschiedlich, da immer die maximal mogliche 
Anzahl von Daten benutzt wird. 
2 . Denormierung von b jk . 

Die Ergebnisse der Identif ikation sind in Datenfiles 
gespeichert mit folgenden Inhalten 

fcjjj.mat - enthalt folgende Parameter, die aus der 

Modellidentif ikation ermittelt wurden: 



V - nicht normierte Clusterzentren 

v - normierte Clusterzentren 

alpha - Locality parameter (hier alpha=0) 

d - Clusterrichtungsparameter 

fuzy - Parameter (hier fuzzy=l) 

u - Zugehorigkeitswerte der Daten zu den 

Clustern nach der Identif ikation 
m - Fuzzyness parameter 

c - Anzahl der Cluster 



Ijjj - Konstante Bias-Werte des Modells 

- Konstante Bias-Werte des Basismodells 
jjj - Jacobians fur alle Cluster in einer Matrix 

zusammengef aSt 

jacob.mat - Jacobians J1...J10 fur jedes Cluster getrennt 
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jacob__b.mat - Jacobians J1...J10 fur jedes Cluster getrennt 

fur das Basismodell 

Die Bettgestalt wird mit Hilfe dreier Parameter, und zwar der 
5 Breite (width) , der Hohe (height) und der Exzentrizitat 

(eccentricity) wiedergegeben . Durch den Operator erfolgt eine 
Zuordnung der von der Kamera wiedergegebenen Bettgestalt zu 
einem ausgewahlten Bettmuster. Nach der optischen Auswahl 
werden intern die drei Bettparameter berechnet und abgespei- 

10 chert. Diese Daten werden wie gewohnliche ProzeSoutputs 

behandelt und werden sowohl zur Prozefiidentif ikation als auch 
zur Optimierung gleichberechtigt mit den anderen Prozefi- 
outputs herangezogen. Allerdings werden diese Daten seltener 
(ca. 60 bis 120 min) als die anderen ProzeSoutputs (miniit- 

15 lich) bestimmt. 

Die Identif ikation der Bettgestalt wird vom Operator durch- 
gefuhrt. Dem Operator wird eine Anzahl moglicher Bettgestal- 
ten als Muster angeboten. Diese Muster erhalten Koordinaten, 

2 0 anhand derer eine zahlenmaSige Benennung eines Musters aber 
auch die zahlenmaSige Angabe eines approximierten Zwischen- 
wertes moglich ist . Die vorliegenden Muster vergleicht der 
Operator mit der Bettgestalt, die er anhand eines Kamera - 
bildes mit Hilfe einer Kamera identif iziert . Hat der Operator 

25 eine Ahnlichkeit zwischen dem Kamerabild und einem Muster 
gefunden, so gibt er in das Advisorysystem die Koordinaten 
desjenigen Musters ein, von dem er annimmt, daS es am 
ahnlichsten mit der im Kamerabild dargestellten angezeigten 
Bettgestalt ist. 
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Ein Kamerabild mit einer erkennbaren Bettgestalt 4 0 ist in 
Figur 4 wiedergegeben . Figur 5 zeigt Beispiele verschiedener 
Bettgestaltmuster . Anhand Figur 5 wird ein Programm (opera .m 
zur Mensch-Maschine-Kommunikation mit Bedienoberf lache in 
WINDOWS -Format fur einen Recovery-Boiler gemafi Figur 1 unter 
Verwendung des vorstehend erlauterten Vorgehens verdeutlicht 

Das Programm opera. m realisiert ein User- Interface zur 
Bestimmung der Bettparameter Hohe, Breite und Exzentrizitat, 
d.h. der (Unsymmetrie . Dem Operator werden verschiedene Bett- 
formen auf dem Feld 61 der Bedienoberf lache 60 in Figur 5 
angeboten, die dieser dann mit der Maus auswahlen und als 
Datenelemente in den Datenvektor Y__FILT eingeben kann. Der 
Operator fungiert somit als "intelligenter Sensor". Sein 
individueller Eindruck von der Bettform wird aufgrund des 
Kamerabildes 40 aus Figur 4 als Datentripel BED_WIDTH / 
BED_HEIGHT, BED_EXC abgespeichert und bei der Optimierung 
beriicksichtigt . 

In opera. m werden mdgliche Bettformen generiert, die der 
Operator individuell auswahlen kann. Die Grundformel zur 
Berechnung der Bettform ist 

y=h*exp (-abs (x/w) . *n) (14) 

mit 

w - Weitenparameter 
h - Hohenparameter 
x - x-Koordinate 
y - y-Koordinate 
n - Exponent 
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Es wird eine symmetrische Kurve erzeugt . Zur Erzeugung von 
asymmetrischen Bettformen wird die Formel gemaS Gl . (14) in 
einen linken und einen rechten Teil geteilt , so daS mit 
einem zusatzlichen Exzentrizitatsparameter auch eine asymme- 
trische Bettform produziert werden kann. Mit vertikalen und 
horizontalen Sliders am Fenster 61 in Figur 5 lassen sich 
Weite, Hohe und Exzentrizitat stufenlos einstellen. Ein 
zusatzliches Fenster 62 zeigt die Parameter der gewahlten 
Bettform an. Mit einem Apply-Button 63 wird die gewahlte 
Bettform als gultig liberaommen. Mit dem Setpoint -Button 64 
kann eine gewunschte Bettform eingegeben werden. Mit dem 
Close-Button 65 wird das Fenster geschlossen. 

Anschliefiend wird die Optimierung des Reglerentwurf s (control 
design) behandelt, wobei ein geeignetes Prozefimodell voraus- 
gesetzt wird: 

a) ProzeSmodell : 

Die Optimierung geht wiederum entsprechend der bereits 
abgehandelten Identif ikation mit Gl . (7) von einem nicht- 
linearen, aber statischen Zusammenhang zwischen Inputs wund 
Outputs y des Prozesses aus : 

^ = /„/(")/ (6') 
deren Approximation bereits erlautert wurde . 

Im Programm wird mit MU(i) bezeichnet und die Vektorenj; 

and u werden durch ihre gefilterten Werte y filt und w^ersetzt. 
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b) Optimierung 

Die Optimierung benutzt die folgende zu minimierende 
Gutef unktion 



J = iy fi n - ysetp) T Q(yjm - y^ ) + - ) T - ) + */i min + -/i max (15) 

mit 

J - Kostenfunktion 

y fiu - n x 1 gefilterter Ausgangsvektor 

y«t P - n x 1 ausgangsseitiger setpoint Vektor 

u fldv - m x 1 Eingangs Advicevektor 

u fih - m x 1 gefilterter Eingangsvektor 

Q - n x n gewichtete Matrix fur Ausgang 

R - m x m gewichtete Matrix fiir Einga ige 

Jl min - Kostenfunktion fiir die kleinsten Beschra nkunge 

Jl max - Kostenfunktion fiir die groflten Beschrankungen 

Es wird angenommen, dafi m < n. Das Optimum wird ermittelt 
durch Gradientenabstieg. Fur jedes Clusterzentrum wird ein 

Vektor Aw'^berechnet , der die Gutef unktion J in Richtung 
ihres Optimums f uhrt : 



^ m -K." m - K . ( ^.^ + *Ly L (16) 

du fiU Wfih «V, du fil, 



Hieraus ergibt sich 

Ak'-* = -K ■ (J e ' T ■ Q • 0> - y tap ) + R ■ ( Uadv - u fiU ) + J e ' T -Y w ) ( 17 ) 

mit * = (1 («)[•// • 2 - (17a) 

1 = 1 

X bestimmt die Lemrate bei der Optimierung. 
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ist ein Vektor, dessen Elemente aus den Output Soft 
Constraints berechnet werden und die nichtlinear von y 
abhangen. Im zugehorigen Prograrnm selbst wurde die Q-Matrix 
durch Ql=abs (sgn (Q) ) ersetzt, urn eine doppelte Gewichtung, 
5 sowohl im Setpoint als auch in der Kostenf unktion zu ver- 
hindern . 

Die Gesamtanderung fur « aA ,ergibt sich zu 

c 

&U adv = -Ai£* • (18) 

1=1 

Anhand des dif f erentiellen ProzeSmodells 

10 

c 

y moi e, (* + !) = y mod e, (*) + Z **t («) • Jj • A "o* (19) 



15 



wird ein u adv (k + 1) = u adv (k) + Au adv (k + 1) solange berechnet, bis die 
Anderungen eine bestimmte Grenze £ unterschritten haben 

y mode/ 

(k)\<s (20) 

beziehungsweise eine bestimmte Anzahl von Iterationen 
erreicht ist. 

Da jedes Outputelement von y nicht unabhangig von den 
Inputelementen von u gesteuert werden kann, ist nicht zu 
erwarten, dafi bei Erreichen des globalen Optimums 

25 J = min auch samtliche Setpoints in y setp erreicht werden. Es 
ergibt sich daher immer ein KompromiS zwichen den zu er- 
reichenden Teilzielen. 
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Auf der Basis der in der Optimierung berechneten Inputs U_ADV 
und dem gerade aktuellen gefilterten Prozefiinput U_FILT sowie 
dem aktuellen gefilterten ProzeSoutput Y__FILT wird uber das 
ProzeSmodell ein Wert Y_EST vorhergesagt, der zusatzliche 
Inf ormationen liber Beschrankungen aufierhalb der Optimierungs- 
schleife berucksichtigt . Y_EST ist derjenige ProzeSoutput, 
der sich ergeben wurde, wenn man unter den obengenannnten 
Beschrankungen den Advice U_ADV implementieren wurde. 

Insgesamt ergibt sich also eine geeignete Operatorberatung 
(Operator- Advice) : Ein Operator -Advice wird dann gegeben, 
wenn 

1. die Stationaritat des Prozesses gestort ist, 
beispielweise durch Parameterf luktuation oder sprunghaften 
bzw. drif tf ormigen Storungen an den ProzeSeingangen, und 
diese durch Veranderung der ProzeSeingange wiederher- 
gestellt werden soli, 

2. ein neuer Sollwert oder eine neue Gewichtung der Ausgange 
gefordert wird und der ProzeS von einem Arbeitspunkt in 
einen neuen Arbeitspunkt gesteuert werden soli. 

Beide Punkte haben Auswirkungen auf die gef orderte Anderung 
des StellgrdSenvektors u . Da es sich nicht urn eine auto- 
matische Regelung handelt, und da der Operator nicht durch 
standige Advices uberlastet werden darf , werden Advices nur 
ab einer bestimmten Schwelle der StellgoSenanderungen 
gegeben. Diese Schwellwerte werden prozeSangepaSt festgelegt. 



Ein sog. Pattern Learning wurde implement iert , urn das stati- 
sche Prozefimodell, das aus einer grofien Menge von Daten ge- 
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wonnen wurde, entsprechend dem aktuellen Prozefiverhalten zu 
adapt ieren. Hierbei wird angenommen, daS die einmal identi- 
fizierten Input-Clusterzentren weiterhin bestehen bleiben. 
Das bedeutet, da£ die Clusterzentren nicht adaptiert werden. 
Zu jedem Zeitpunkt werden gleitende Mittelwerte 



DELTA 7=72-71 = 



1 



N JAC 



MAXH1ST-X 



MAXHIST-2 



2 r_cuR(i)- 5Xc«R(o 

i=MAXHlST-N _JAC <=MAXH1ST-N JAC-} 



(21) 



10 DELTA U=U2-U1 = - 



1 



N JAC 



MAXH1ST-I 



MAXHIST-2 



]T U_ CUR(i) - ]T U_ CUR(i) 



i - MAXHJST- N JAC 



i = MAXH1ST- N _JA C- \ 



(22) 



berechnet. Dabei sind U_CUR und Y_CUR ungefilterte Ein- bzw. 
Ausgange, N_JAC die Anzahl der zu mittelnden Werte und 
MAXHIST eine Zahlvariable . Danach werden aus dem Modell die 
15 Werte 



delta _y = ]T MU(i) • J* (0 * DELTA _ U 



(23) 



20 



berechnet . 



Das Adaptionsgesetz fur eine neue Jacobian J NEW lautet : 



J NEW i = J< c + LearRat • MU(i) ■ ^LTA_Y - delta _y) • DELTA_U T 
~ DELTA_ U T ■ DELTA_ U + Spoot 



25 mit der Adaptionsrate LearRat und der Konstante Spoot. Mit 
der Einfugung des Multiplikator MU(i) werden die Cluster- 
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zentren entsprechend ihrer Anregung adaptiert, wodurch die 
sog. "Persistence of Exitation" berucksichtigt wird. Um ein 
zu grofies Abdriften der neu berechneten Jacobian von der 
urspriinglich berechneten Jacobian zu verhindern, wird die 
Norm der Differenz DELiTA_J_i = J_NEW__i - J_B_i berechnet, die 
nicht grofier sein darf, als eine vorgegebene Schranke THRES. 
Wird diese Schranke erreicht, so andert sich die Jacobian in 
diesem Schritt nicht . 

Schliefilich ist auch ein sog. Batch Learning moglich: 
Um Arbeit spunktveranderungen wirksam abfangen zu konnen, wird 
zyklisch eine Neuberechnung des Prozessmodells durchgef uhrt . 
Hierzu werden die Daten, aus denen das urspriingliche Grund- 
modell berechnet wurde, mit aktuellen Daten angereichert und 
mit diesem erweiterten Datensatz sowohl eine neue Cluster- 
bestimmung als auch eine neue Modellberechnung durchgef uhrt . 

Eine zusatzliche Option auf eine Operatorberatung ist, einen 
Advice dann zu geben, wenn der Absolutwert der geforderte 
StellgroSenanderung einen bestimmten Prozentsatz des Abso- 
lutwertes der StellgroSe angenommen hat . 

Der Advice wird wie f olgt berechnet : 

1. Einlesen der aktuellen ProzeSinputs und -outputs 

2. Filterung der aktuellen Prozefiinputs und -outputs 

3. Bei Bedarf Eingabe der Bettform 

4. Berechnung des Advices aus der Optimierung und einer 
entsprechenden Pradiktion 

5. Test, ob die geplante Adviceanderung unterhalb einer 
vorgeschriebenen Schranke liegt, womit die Empf indlichkeit 
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definiert ist.Dabei ergeben sich zwei Moglichkeiten : 

- Wenn ja — » keine Adviceanderung 

- Wenn nein — > Adviceanderung 

6. Aussetzen des Advices fur eine vorgeschriebene Zeit , wenn 
5 der Operator den Advice implement iert hatte . 
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Patentanspriiche 

1 . Verf ahren zur Prozefif uhrung einer verf ahrenstechnischen 
Anlage, insbesondere eines Ruckgewinnungskessels (sog. 
Recovery-Boiler) bei der Zellstof f herstellung, wobei aus fur 
Eingange und Ausgange anfallenden ProzeSdaten Signale ermit- 
telt und zur Steuerung und Regelung des Prozesses verwendet 
werden, mit folgenden Verf ahrensschritten : 

- Die Rohdaten werden in einer Einheit zur Datenvorver- 
arbeitung aufbereitet , 

- die aufbereiteten Daten werden einer Systemidentif izierung 
unterzogen, 

- voiti Operator werden gewiinschte Zielwerte der Ausgange, 
vorzugsweise ProzeSsollwerte , eingegeben, 

- darnit wird eine Kontroll-/Optimierungsstruktur abgeleitet, 
die eine Veranderung der Eingange vorschlagt, 

- die strukturierten Daten werden in den Regelprozefi ein- 
gegeben . 

2. Verf ahren nach Anspruch 1, dadurch gekenn- 
zeichnet , da£ zur Systemidentif izierung 
datengetriebene Methoden verwendet werden. 

3. Verf ahren nach Anspruch 1, dadurch gekenn- 
zeichnet , dafi zur Systemidentif izierung Fuzzy- 
Methoden, insbesondere eine Fuzzy-Clusterung, angewandt 
werden . 

4. Verf ahren nach Anspruch 1, dadurch gekenn- 
zeichnet , daS zur Systemidentif izierung statisti- 
sche Methoden verwendet werden. 
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5. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet , daS das Prozefimodell aufgrund von aktuel- 
len ProzeSdaten im Sinne einer Adaption nachkorrigiert wird. 

6 . Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet , dafi vom Operator die Gewichtung der zu 
erreichenden Zielwerte einstellbar ist . 

7. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet , daS dem Operator ein Vorschlag fur ein zu 
implement ierendes Steuersignal unterbreitet wird. 

8. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet , daS dem Operator entsprechend dem Steuer- 
vorschlag und dem Prozefimodell eine Vorhersage der zu erwar- 
tenden ProzeSreaktion angeboten wird 

9. Verfahren nach Anspruch 1, wobei ein Recovery-Boiler mit 
einem Schmelzbett (char bed) variabler Gestalt vorhanden ist, 
dadurch gekennzeichnet, daS 

- die Gestalt des Schmelzbettes mit einem Katalog typischer, 
automatisch parametrierbarer Bettgestalten verglichen wird, 

- durch optischen Vergleich des Schmelzbettbildes mit den 
angebotenen Mustern eine Grobklassif ikation der Bettgestalt 
vorgenommen wird, 

- eine Feineinstellung der Parameter des ausgewahlten Musters 
dieses so verandert, daS es dem Schmelzbild am nachsten 
kommt , 

- die Parameter des ausgewahlten und f eineingestellten 
Musters in die Datenbasis automatisch eingegeben werden, 
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- die eingegebenen Bettparameter als aktuelle Bettdaten 
ve r arbe i t e t we r den . 

10. Verfahren nach Anspruch 9, dadurch g e - 
kennzeichnet , daS das Schmelzbett durch eine 
Kamera erfafit und als Kamerabild dem Operator auf einem 
Monitor vermittelt wird. 

11. Verfahren nach Anspruch 9, dadurch g e - 
kennzeichnet , daS der Vergleich und/oder die 
Feineinstellung durch den Operator erf olgt . 

12 . Vorrichtung zur Durchf uhrung des Verf ahrens nach Anspruch 
1 oder einem der Anspruche 2 bis 11, vorzugsweise zur ProzeS- 
f uhrung eines sog. Recovery-Boilers als verf ahrenstechnische 
Anlage, wobei der Anlage wenigstens eine Steuer- und/oder 
Regeleinrichtung zugeordnet ist, dadurch g e - 
kennzeichnet , dag ein Regelkreis mit einer Ein- 
richtung (Advisory- System 30) zur Generierung von Signalen, 
fur die wenigstens eine Steuer- und/oder Regeleinrichtung 
einschliefilich Regelkreis vorhanden ist, mit dem einem 
Operator (OP) die neuen Steuer- und Regelsignale mitgeteilt 
werden . 

13. Vorrichtung nach Anspruch 12, dadurch g e - 
kennzeichnet , daS das Advisory- System (3 0) ein 
Interface (60) beinhaltet, an dem vom Operator (OP) das Bild 
(40) des Schmelzbettes mit den vom Interface (60) angebotenen 
Bettformen verglichen und durch Feineinstellungen an die 
aktuelle Bettform angepafit wird. 
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